Generation of Superposition States and Charge-Qubit Relaxation Probing in a Circuit 
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We demonstrate how a superposition of coherent states can be generated for a microwave field 
inside a coplanar transmission line coupled to a single superconducting charge qubit, with the 
addition of a single classical magnetic pulse for chirping of the qubit transition frequency. We show 
how the qubit dephasing induces decoherence on the field superposition state, and how it can be 
probed by the qubit charge detection. The character of the charge qubit relaxation process itself is 
imprinted in the field state decoherence profile. 
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I. INTRODUCTION 

Generation of non-classical states of the electromag- 
netic field has been one of the most pursued problems in 
quantum optics pQ. This interest was reinforced in the 
last years due to its potential applicability in quantum 
information 2\. Aside from the experimental effort to 
generate single photon states, the generation of super- 
position of coherent states has been for a long time a 
central task, with some successful attempts. Superposi- 
tion of coherent states, commonly known as Schrodinger 
cat states, were firstly generated in a superconducting 
microwave cavity field interacting with flying Rydberg 
atoms and its decoherence time was measured in [3J also 
employing Rydberg atoms sequentially interacting with 
the field following a previous theoretical proposal [4j. 
More recently these states were generated in a propa- 
gating field by photon-subtraction from a Gaussian state 
[5]. Indeed, special superpositions, known as odd and 
even coherent states have well defined parities and can 
have immediate application to encode quantum bits in a 
robust way [5] . In the last few years a new technology for 
coupling superconducting qubits to coplanar waveguides 
has been developed to an outstanding level [7] . Many dif- 
ferent tests have been developed and interesting quantum 
optical experiments have been implemented in supercon- 
ducting circuits [5]. Then, a simple question could be 
posed on how to generate a superposition of coherent 
states in this kind of system, and subsequently how to 
probe its decoherence due to dissipative effects of the 
qubit or the field. To our knowledge the only proposals 
in that direction employ a SQUID-type two-level system 
coupled to a microwave field !) . and a micromechanical 
resonator coupled to a Cooper-Pair Box [TU] (See Ref. 
[9] for a complete list of connected proposals in distinct 
technologies). However none of those proposals bear a 
resemblance with the neat experiment in [3J. 

In this paper we show how to proceed to generate a 
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FIG. 1. (a) Schematic setup, with the central transmission 
line (resonator) capacitively coupled to the source and drain, 
and capacitively coupled to a SQUID-type qubit. (b) Variable 
energy levels of the qubit (0 < n g < 1), with an external 
classical magnetic flux & x (t). The resonator field is always 
blue detuned from that transition. 



superposition of coherent states of a microwave field in a 
transmission line resonator through the interaction with 
a single superconducting charge qubit controlled by an 
external single classical magnetic pulse. In contrast to 
the proposal in Ref. [9], which employs of a switchable 
interaction between the qubit and the microwave field, we 
assume a continuously varying magnetic pulse, and show 
how the generation can be almost deterministic in that 
way. The present proposal has the advantage to resemble 
the experiment in [3J in a realistic ground, allowing thus 
for decoherence probing through sequential qubit mea- 
surements and sequential classical magnetic pulses. We 
show that the qubit states dephasing is the most relevant 
dissipative effect leading to decoherence of preselected 
field states. Employing actual experimental parameters 
we show that this proposal is compatible with present 
technology, and could be actually implemented. 



II. TIME DEPENDENT FIELD-QUBIT 
COUPLING 

Cavity quantum electrodynamics in superconducting 
circuits offer a exquisite playground for quantum infor- 
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mation processing, and has provided the first coherent 
coupling between an "artificial atom", the charge qubit, 
and a field mode of a resonator [7J [TT]. Mappings of 
qubit states [T2], and also tests for fundamental prob- 
lems, such as the Purcell effect [T5] and photon number 
state resolving 13j have also been achieved. The setup 
employed in all those remarkable experiments is shown 
in Fig. |la[ where a niobium transmission line resonator 
is capacitively coupled to a source (on the left) and to a 
drain (on the right). The resonator is also capacitively 
coupled to the charge states of a SQUID [TT] . The advan- 
tage of employing a SQUID is that the charge states can 
be addressed and manipulated in such a way to be set 
close or far from resonance with a given resonator field 
mode by an externally applied classical magnetic flux. By 
considering only the ground and the first excited states 
near the charge degeneracy point, the superconducting 
device can be well approximated by a two level system 
(Fig. lb), here addressed as a qubit. In this regime the 
Hamiltonian |llj describing a quantized electromagnetic 
field mode coupled to the charge qubit is given by 



H = tkd a' a + Hi 



where 



He 



'-B.a. 



-B x a x , 



(1) 



(2) 



is the Hamiltonian for the qubit-ficld joint system, where 
B z = AE c {l-2n g ), being ,E C = e 2 /2{C g +C ] ) = e 2 /2C s , 
the single electron charging energy with C g as the gate 
capacitance ( associated to the accumulated charge in 
the island capacitively coupled to the central transmis- 
sion line resonator) and and Cj the Josephson capaci- 
tance, respectively. n g = C g V g (t)/2e and at the center 
of the transmission line resonator the voltage is given 

by V g (t) = \J~^{a i + a), where L is the length and 
c is the capacitance density of the transmission line. 
Here, a (aJ) is the usual annihilation (creation) opera- 
tor for the resonator field second mode of frequency id, 
a x = |0) (1| + |1) (0| and a x = |0) (0| - |1) (1|, involv- 
ing the ground, |0), and first excited, |1), charge states 

of the superconducting device. B x = 2Ej cos n 

is the Cooper-pair tunneling energy for two Josephson 
junctions, $o = hc/2e is the quantum of magnetic flux, 
$ x (i) is a time dependent classical magnetic flux exter- 
nally applied to the SQUID. All other quantities are typ- 
ical constants characteristic of the device [TT] . Collecting 
all terms the Hamiltonian, Eq. reads 



H = Huxr a + Ej cos 



$»(*) 
<J>n 



+ Hga z (a t + a) (3) 



being g = (eC ' g / Cs) \J Hut / Lc the coupling between the 
"artificial atom" (qubit) and the resonator mode. In Eq. 
^ we neglected, as usual, the term — (e 2 /Cs)er 2 corre- 
sponding to a DC voltage shift. 

When the atomic Hamiltonian is diagonalized, the 
first two energy levels as a function of the gate charge 



CgV g /2e are described in fig. lb where the verti- 



cal axis represents the energy and the horizontal repre- 
sents the gate charge which is limited by the gate volt- 
age. Changing the basis through a rotation, a z — > a x 
and a x — > —&z, and going to the rotating frame with 
the field frequency, uj, through Rf = exp \icdt(a z + at a)] , 



gives 



d 
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with the transformed Hamiltonian given by 



dt 



1> 



H = cr. Ej cos 



— froj^J + hgd x (aJ + a) (4) 



where a x = er + exp(2iwt) + cr_ exp(— 2iut), a = 
aexp(-iujt) and a' = (V exp(iwt), with a x = \— ) (+| + 
|+) <-| and a z = |-) (-| - |+) (+|. This new Hamil- 
tonian Q, is analogous to the usual one for atom-field 
dipole interaction. The appropriate experimental regime 



is hg « (^Ej cos 



.Ml 



hu; J « A, where A is the 
superconducting energy gap. In the following we shall 
consider a specific time dependent external classical mag- 
netic pulse $ x (t) applied to the SQUID, in order to bring 
the two lower states closer to resonance with the res- 
onator field. As we will show, during the pulse, the field 
accumulates a phase conditioned to the qubit state. For 
that we consider the resonator field as being always far 
blue-detuned from these two lower atomic states, and so 
we must keep the counter-rotating terms 13J . With that 
in mind we derive the formal solution for the evolution 
operator 
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xV I {t 1 )V I {t 2 )...V I {t n ), 
with V '(i) ee Ul {ga x (a f +a)}U , and 



dt. 



(5) 



Uq = exp 



dt 



Ej 



For k B T « Ej « E c « A, a perturbative ap- 
proach up to second order in g in Eq. ^ is sufficient 
to describe the dynamics of the system. For that the 
temperature must be as low as 30mK, which is con- 
sistent with this kind of experiment. In what follows 
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we consider the solution of ([5]) through time-dependent 
perturbation with Hq — a z ( Ejcos 



Huj ] and 



Vf = hg (a + + a ) (fit + a), with H » Vf , since we 
are interested in the radiation field blue detuned from the 
qubit transition, which on its turn oscillates with time. 
For our numerical calculations, we use the following clas- 
sical magnetic flux 



$»(t) 



cos [at + tp] , 



(6) 



where A — 0.7 is a strength constant, and a = 8ttx 10 6 Hz 
is the frequency of the classic magnetic pulse applied 
to the qubit in accordance with experimental values. 
By further assuming typical experimental parameters to 
reach k B T « Ej « E c « A, for T « 30mK, we 
set k B T « 3/ieV, Ej/h w 15,9 x 10 10 Hz, E c = 250/zeV 
and A « 458, 3/ieV. The frequency of the field in the 
resonator, lu = 90 x 10 10 Hz, is also compatible with 
the experiments, for which the resonator quality factor 



Q = 



10 4 - 10 6 is achievable [IS]. 
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FIG. 2. Phase of the resonator coherent field, due to the 
time dependent interaction with the qubit prepared in a) the 
ground state, and b) the excited state. The pulse oscil- 
lates for a half period (At ~ 7.5ns) with frequency 8irx 10 6 Hz. 



following approximation for short time, 1 + 9±(t)a'a 
e [e±{tWa] ( gee the Appendix) so that 



[1 + e±(t)<fa] |±) \a) -> e -5M 2 /±(t) | ± ) 



9±(*) 



(7) 



where /±(t) = [1 - e 2tfe(<?± (*))], ± ( t ) = i + ; being 



' dt e™^ 1 dt ' n (*')) = F™*i 

±2i(/ * 2 dt"n(t"))±itut 2 ^ 
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G±(t) 



dt ^{So 1 dt n{t'j)±iwtx 
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x / di 2 e 
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with fl(t) = uj it cot 

order terms from Eq. 



as given by the second 



From Eq. Q it is possi- 
ble to understand that as the qubit is brought closer to 
resonance with the resonator field, it will imprint an ac- 
cumulated phase on it, given by lm[0±(t)] conditioned on 
the qubit state |±). In Fig. [2]we depict the numerical re- 
sults for those two conditioned accumulated phases. We 
see that practically only when the qubit is in the state 
|— ) the phase in a is changed. With an appropriate ac- 
cumulation of — 37T, as shown in Fig. [2^,, it is possible to 
create a state |— a) if the qubit is initially in the |— ) state. 
We have consistently checked that this approximation is 
indeed very good, not only for small a, if we respect 
a balance between the field intensity and the operation 
time. Moreover, we observed that around the time of op- 
timal phase accumulation, t op = 7.5ns, the real part of 
9±(t), related to damping or amplification, is negligible, 
(ra 10~ 3 - 10~ 4 ), so that in Eq. @, f ± (t) = at the 
instant of measurement as shown in the Appendix, and 
shall not be considered from now on. 

To generate the superposition of coherent states for 
the field, we note that due to the low temperature of 
the system it is easy to prepare the qubit state initially 
in |0). Obviously the decoupled coherent state of the 
resonator and the qubit state may be written as |0) ® 
| a) = (|— ) + |+)) ® \a)/V2. If the aforementioned pulse 
is applied to the qubit, coupling it to the resonator field 
through the evolution given by Eq. (j5j|, we shall have 



III. ACCUMULATED PHASE AND 
CONDITIONAL GENERATION OF 
SUPERPOSITION STATES 

In order to depict our results in a more convenient way 
we choose that the field in the central line resonator is 
prepared in a coherent state \a) with an average number 
of photons smaller than one, so that we can make the 



U{t +At,t )\0)®\a)=> !-t '— — , 

V 2 

Pulse 

(10) 

or [|0) <g> (\-a) + \a)) + |1) <8> (|-a) - |a))]/2, which is 
an entangled state between the qubit and the resonator 
field. Consequently, the resonator field can be left in an 
odd or even superposition of coherent states depending 
on the detection of the qubit state with a single electron 
transistor H4l. 
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IV. FIELD STATE DECOHERENCE AND 
QUBIT RELAXATION PROBING 

Certainly, the exact preparation of such a state is com- 
promised by external noise. In contrast to experiments 
with microwave fields and Rydberg atoms, dissipative ef- 
fects are most noticeable for the qubit states, which can 
flip from the ground state to the excited one and vice 
versa due to thermal effects and inductive coupling of 
the qubit to the external circuit. While the relaxation 
time of the qubits are of the order of 10 -6 s , the relax- 
ation of the resonator field is of the order of 10~ 3 s and 
can, in principle, be neglected. If compared with the 
time of the pulses for the accumulated phase of the field, 
At w 7.5 x 10~ 9 s , the relaxation of the atom is negligi- 
ble as well, but certainly will be important if any further 
manipulation is to be executed. So in fact the effects of 
dissipation are more relevant after the pulse is applied, 
i.e., when $ x (t) = 0, and will affect the probability to de- 
tect a given qubit state, and consequently the generation 
of an appropriate field. To understand the effects of noise 
in the system, we couple the qubit two-state to a bath 
of harmonic oscillators in an adaptation of the standard 
spin-boson model with Ohmic dissipation []3J [T71 [T5] . 
through the Hamiltonian H = Hs + Hab + Hr- Here H$ 
is given by Eq. presented in reference [14], with the 
convention at the degeneracy point |+) = (|0) — |1})/V2, 
and |-) = (|0) + \l))/y/2. Eq. Q can be conveniently 
re writ cn as 



AE 1 1 

H S = - — {7=2 * Z + 7=2 * X ), 



(11) 



where AE = \AB| + B%, and now since <fr x (t) = 0, 



B z = 4 £ c (l 



16^(1 



2%), and B x = 2Ej. Thus AE = 
and conveniently choosing 



2n ff ) 2 , 



0.5 we end up with AE — 2Ej. The coupling of 



the qubit to the bath is given by 



Har = &z ^ A a x Q , 



and 



E 



2„2 
a^a x a 



Pa . lll„X-.r 

2m a 



(12) 



(13) 



corresponds to the bath free Hamiltonian. 

The master equation for the evolution of the two 
level system coupled to a bath in a thermal state 
Pn = cxp [H^/ksT] was extensively studied in the past 
[T71ll8j . We employ solutions corresponding to an Ohmic 
bath [TH [T7J HHj. By taking the previous initial state, 
1 0)| a), the main consequence of the qubit relaxation will 
be on the probability to generate the superposition states 
|0) L = (\a) + \-a))/V2, or = (\-a) - \a))/y/2 af- 
ter the pulse is applied. To turn our description clearer 
we consider the optimal time for phase accumulation as 
our time origin. Remark that since the central trans- 
mission line resonator field is now far detuned from the 



qubit, they remain uncoupled. The only possible correla- 
tion between the qubit and the field is due to their past 
interaction. Taking the result for the density operator 
evolution from Refs. Q31 [T71 [T5] we obtain the following 
joint qubit-ficld state 

P (t) = P (t)\0)(0\ ® |0)(0|x + Pi(t)\l)(l\ ® |1)(1|l (14) 
+Pr(i)|l>(0| ® |l)(0|i + i?(i)|0>(l| ® |0)(1|l, 
where P(o)(t) = Tr{\ (°)) ((°) \p(t)} is the probability to 
detect the qubit in the state Id)), given by 

P(oj(t) = - {tanh(A) + [1 - tanh(A)] exp(-t/r r ) 

± cos(AEt/K) exp(~t/T v )} , (15) 

and Pr{t) = Tr{\0)(l\p(t)} is the transition amplitude, 
explicitly given by 



P T (t) = -ism(AEt/h)exp(-t/T v 



(16) 



where AE = 2Ej, A = AE/2k B T, r r = 
[7r/3Ai?coth(A)//i] _1 is the relaxation time and t v = 
[t j T 1 /2 + 2TTf]kBT/H]~ 1 is the dephasing time, with j3 
0.001 a dimensionless parameter reflecting the strength 
of dissipation [13]. As AE > k B T, tanh(A) w 1, and so 

Po(t) = i[l + cos(AEt//i)exp(-Vv)], (I 7 ) 
Pi(t) = i[l - cos{AEt/h)cxp(-t/T v )}, (18) 

reflecting the probability to detect the qubit in the state 
|0) or |1), respectively, at an instant t after the classical 
magnetic pulse. The entangled state from Eq. ( 14 1 allows 



a probabilistic generation of the field superposition state, 
as well as an indirect probe of the qubit state by the field 
state measurement. 

As can be readily checked, if the qubit detection is 
made right after the pulse is applied, there is a high prob- 
ability that the state |0) L will be generated, in this case 
thus almost deterministically. If the detection is made 
after a long delay time, t sa 3r v , there is 50 % chance 
that the qubit can be detected in |0) and 50 % in |1) 
and thus that the states |0) L or \l) L , be generated, re- 
spectively. Indeed, if a second classical magnetic pulse 
is applied, after a qubit detection, the state generated 
is exactly equal to (14) and a sequential detection after 



a time interval t, allows the inference of those probabili- 
ties directly from the measurement, since the conditioned 
probability for two consecutive detections in the state |0) 
is Po,o(t) — Po(t), and the probability for detecting |0) 
and then |1) is Po i(i) = Pi(t). This contrasts with the 
decoherence time detection in Ref. [3]. Those probabili- 
ties are depicted in Fig. 3, for T = 30 mK. There we em- 
ploy an artificial AE, which is 10 3 times smaller than the 
real one to depict the oscillatory profile. The probability 
to sequentially detect |0) decreases with time, as a signa- 
ture of the dephasing of the qubit states and thus on the 
probability to generate \0}l- Then of course the proba- 
bility to generate \0) l increases with time. The decaying 
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and charge qubits. One immediate example is the Bell- 
like state Eq.(8), for t « r,,, which may directly allow the 
implementation of quantum communication protocols, or 
indirectly through a possible coupling to other solid state 
qubits (TBJ- We leave this point for future discussion. 
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FIG. 3. Probabilities of charge qubit detections and conse- 
quently of postselected |0}l, or 1 1) z, field states. The deco- 
herence of the preselected field state is given by 2Po(t) — 1. 



oscillatory profile depicted in Fig. 3, which represents 
exactly the dephasing of the qubit states, will surely be 
imprinted in the field states even in the absence of detec- 
tions, and is a signature of the Ohmic character of the 
reservoir. The preselected field state derived from Eq. 



APPENDIX 



(14) is 



pf(t) = t777t{I«)H + I -«><-«! 



N(t) 



(19) 



+ cos(-^-)e~ t/T »° [\a)(-a\ + | - a)(a| 



with N(t) = 2[l + cos(AEt/h)e~ t / T ^e- 2 ^ 2 }. So, the de- 
phasing causes the decoherence of the state pf(t), whose 
behavior is exactly the one depicted in fig. 3, and can be 
inferred as being 2P (t) — 1. We remark that in the same 
sense the decoherence of the field can be used to probe 
different noise characters on the qubit. Here we analyzed 
the effect of an Ohmic bath. Would that be a super- 
Ohmic or sub-Ohmic its signature would be imprinted in 
the field decoherence. 



V. DISCUSSION 

We have discussed the possibility to generate superpo- 
sition states of the field due to a classical magnetic pulse 
which causes a chirping of the frequency of the qubit, 
bringing it closer to resonance with the resonator field 
mode. At low temperatures the dissipative effects of the 
field mode are almost negligible during the time inter- 
val the superposition is generated, while the dephasing 
of the qubit becomes the relevant source for decoherence 
of the state of the radiation field. In this case the qubit 
detection probability affects the probability for genera- 
tion of the field superposition states |0)l or \1)l- Those 
states are quite robust against dissipation |6j, and there- 
fore may be quite relevant for quantum information pro- 
cessing involving a hybrid form of qubits involving field 



In this appendix we show how Eqs. ([7|-([9]) are derived. 
For that we have to perturbatively compute how the evo- 
lution operator, U(t,to), shown in Eq. ^ acts over an 
initial state |0) |a), 



U(t,t o )\0)\a)=U{t,t o ) 



V2 



(A-l) 



Due to the parity of the classical flux $> x (t) applied over 
the qubit, Eq. Q, for a "fast" pulse, i.e., a pulse with 
duration equal to or shorter than half period of oscilla- 
tions (for to — 0), the first order terms in Eq. ^ are 
negligible and we consider only terms of second order. 
So that, 

U(t,t ) |+) \a) « -g 2 [a^aF + (t)+aa^G + (t)] |+) \a) , 

U(t, t ) |-) \a) w -g 2 [aa}G_ (t) + a^aF_ (t)] |-) \a) , 

where F± (t) and G± (t) are given by Eq. Q and Eq. (j^, 
respectively, 

f dhe™^ 1 dt ' f2 (*')) = F™*i 
lo 



F±(t) 



dtoe 



±2i(/ * 2 dt" fi(i"))±uot 2 



G± (t)= /' t ^ ie ^(/ tl *'o( t '))±™t 1 



±2 



i(/ * 2 dt ' U(t")J^fiwt 2 



corresponding to terms of second order of the operator 
U(t,t ) from Eq. ((Hj). 

The second member of the evolution equations above 
can be respectively rewritten as 

~g 2 {G+(t) + [F + (t) + G+(t)} a^a} \+) \a), 

-g 2 {G_(t) + [F_(t) + at a } I-) \a), 

or alternatively as 

, '- [ '^a^\+)\a), (A-2) 



-g 2 G+(t)^l + 
-g 2 G-(t)(l + 



G + {t)_ 
F-{t) 
G-(t) 



i*a\\-)\a). (A-3) 
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FIG. 4. The function f±(t) operating on the states |±). 



The terms —g 2 G±(t) give only contributions to a global 
phase of no implication and without any loss of generality 
are being neglected. For short time and small average 



number of photons, 



1 



F±(t) 
G±(t) 



a) a is kept small enough 



so that the following approximation can be employed 



F±(t) 



F±(t) 



G±(t) 



and thus 

[1 + 6 ± {t)a)a] |±) | a) -> e^M'/M |±) 



a) a 



J±(t) 



1 



in Eq. 

G±(t)' 



7b where f±(t) = [1 



o 2Re(0±(t))l 



=(*) 



This approximation has to be taken with cau- 
tion. As shown in Fig. Ek.) , /_ (t) increases to large (neg- 
ative) values, achieving its maximal value around 4.5ns 
meaning that at those times the approximation in Eq. ^ 
is not good enough. However closer to the time of opti- 
mal phase accumulation, t op = 7.5ns, in Fig. 
decreases rapidly to (w 10~ 3 — 10~ 4 ), meaning that the 
accumulated real part of the pulse is negligible well be- 
fore the end of the pulse. On the other hand the function 
f+(t) = 1 — e 2Ro ^ e+ ^\ shown in Fig. (J^), is always very 
small, of order 10 -3 — 10~ 4 , during the time of the pulse. 
Thus showing the validity of the approximation in Eq. 
Q, is consistent for the time of the pulse. 

Now at the optimal time, t op , since f±(t op ) ~ 0, Eq. 
([7]) can be effectively written as 



l±> 



J±(ta P ] 



(A-4) 
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